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ABSTRACT 

A proposed method for dealing with foreground emission in upcoming 21-cm observations 
from the epoch of reionization is to limit observations to an uncontaminated window in 
Fourier space. Foreground emission can be avoided in this way, since it is limited to a wedge- 
shaped region in A:||, k^_ space. However, the power spectrum is anisotropic owing to redshift- 
space distortions from peculiar velocities. Consequently, the 21-cm power spectrum measured 
in the foreground avoidance window—which samples only a limited range of angles close to 
the line-of-sight direction—differs from the full redshift-space spherically-averaged power 
spectrum which requires an average over all angles. In this paper, we calculate the magnitude 
of this “wedge bias” for the first time. We find that the bias amplifies the difference between 
the real-space and redshift-space power spectra. The bias is strongest at high redshifts, where 
measurements using foreground avoidance will over-estimate the redshift-space power spec¬ 
trum by around 100 per cent, possibly obscuring the distinctive rise and fall signature that is 
anticipated for the spherically-averaged 21-cm power spectrum. In the later stages of reion¬ 
ization, the bias becomes negative, and smaller in magnitude (< 20 per cent). 
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1 INTRODUCTION 

A major obstacle facing any upcoming measurement of the 21-cm 
signal from the epoch of reionization (EoR) is foreground emission 
from galactic and extra-galactic sources (e.g. Jelic et al. 2010). Two 
separate approaches for dealing with this problem have been pro¬ 
posed: foreground removal and foreground avoidance. Foreground 
removal involves modeling the foreground emission in order to sub¬ 
tract it from the total signal (e.g. Chapman et al. 2013, 2012). In 
foreground avoidance, the idea is instead to use only the parts of 
the signal in k space (Fourier space) that are unaffected by fore¬ 
grounds. 

It has been shown that foreground contamination of the FoR 
21-cm signal will he localized to a wedge-shaped region in cylin¬ 
drical fc||, k± space, at low values of p = A:|j /|k| (Datta et al. 2010; 
Trott et al. 2012; Dillon et al. 2014; Pober et al. 2014). This leaves 
an “FoR window” at high p, where the pristine signal can be ob¬ 
served. Limiting observations to this window has the advantage of 
not requiring any knowledge of the properties of the foreground 
emission. The main downside is that a (potentially large) part of 
the signal has to he thrown away, resulting in lower signal-to-noise. 
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A second downside is that direction dependent quantities, such as 
the /r-dependent power spectrum become hard to measure (Pober 
2015). 

Fven for measurements of the spherically-averaged power 
spectrum (which is not dependent on pt), one must he careful when 
comparing observations to theory. Observations can only be done 
in redshift space, where the apparent positions of emitters are dis¬ 
torted by their peculiar velocities along the line-of-sight, and it is 
well-known that the redshift-space power spectrum, As(k), can 
differ substantially from the real-space power spectrum, Ar(k), 
even after spherical averaging (Mao et al. 2012)^. However, since 
the redshift-space power spectrum is not spherically symmetric, 
limiting a measurement to an FoR window—i.e. including only 
some values of pt —will introduce an additional bias. We then 
have to distinguish between three different power spectra: Aj (fc), 
which is the fundamental one, and the raw output from most sim¬ 
ulations; Agjfc), which is the quantity seen by an observer; and 
Ag window (^)> which is a biased estimator of Ag (k). 


^ We show our results in terms of the “dimensionless” 21-cm power spec¬ 
trum A‘^{k) = P{k )/We use subscript r to denote real-space 
quantities, and s to denote redshift-space quantities. 
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The bias that makes Ag (k) different from Ag (k) has 

been hinted at in the literature (e.g. Liu et al. 2014), but to our 
knowledge, it has never been quantified. Here, we calculate the 
magnitude of the bias and discuss how to correct for it when in¬ 
terpreting observations and simulations. 


2 THEORY 
2.1 The wedge 

The foregrounds that contaminate EoR 21-cm measurements are 
expected to be spectrally smooth, which means they should only 
affect the lowest fey modes. However, because of the frequency- 
dependence of an interferometer’s response, the contamination will 
leak into higher k modes. This effect, known as “mode mixing”, is 
strongest for long baselines (high values of A:^) which have higher 
fringe rates. Mode mixing causes the foregrounds to become con¬ 
fined to a wedge-shaped region in fc||, k± space (Datta et al. 2010; 
Trott et al. 2012; Pober et al. 2014; Dillon et al. 2014; Liu et al. 
2014; Pober 2015). 

The contamination by flat-spectrum foreground sources can be 
shown to extend no farther than the line defined by the following 
relation: 

, _f ■ Q Dm{z)E{z)\ _ ,,, 

V Dh{1 + z)J 

where 9foV is the angular radius of the field-of-view, Du = c/Hq, 
E{z) = ^Hm(l + + Ha and Dm{z) is the transverse co¬ 

moving distance, which for a flat cosmology is equivalent to the co¬ 
moving distance: Dm{z) = Du Jq dz'/E( z'), (e.g. Parsons et al. 
2012; Thyagarajan et al. 2013; Dillon et al. 2014). 

From the slope of the wedge defined in Equation (1) we get 
the minimum value of that can be observed without foreground 
contamination: 

^ ' 

Unfortunately, there is some lack of consensus as to what value for 
0Fov should be used when calculating the wedge slope. The most 
pessimistic assumption would be to include contamination from 
sources on the horizon, and set Spov = 90 degrees, although obser¬ 
vations seem to indicate that some “supra-horizon” contamination 
may occur beyond this limit, possibly since real foregrounds are in 
fact not spectrally flat (Pober et al. 2013). On the other hand, it has 
been argued that it is possible to avoid contamination from sources 
outside the primary beam of the telescope, which would make the 
EoR window significantly larger (Pober et al. 2014). 

In Figure 1, we show /imin as a function of redshift for a few 
different field-of-view sizes. Our aim with this paper is not to pro¬ 
vide a realistic model for the EoR window, so we will simply show 
our results for the optimistic case of /Tmin = 0.5 and the more 
pessimistic case of /imin = 0.95. 

In addition to the wedge, the EoR window is also restricted by 
foregrounds on all low fey modes. This edge is commonly assumed 
to be around k\i = 0.05 Mpc“^ (Dillon et al. 2014). The lowest k 
mode we consider in this paper is A: = 0.09 Mpc“^ which is well 
within the window. First-generation telescopes such as LOFAR and 
MWA are expected to be most sensitive to the EoR 21-cm power 
spectrum on scales around k ~ 0.1 Mpc“^ (Harker et al. 2010; 
Beardsley et al. 2013). 



Figure 1. The smallest fi that is observable using foreground avoidance, as 
a function of redshift for different field-of-view sizes. A 90 degree field-of- 
view corresponds to the horizon limit. 



Figure 2. A slice from the full three-dimensional 21-cm power spectrum 
at 2 = 10, smoothed with a Gaussian kernel for visual clarity. The shaded 
black area shows the region where /i < 0.5. The white region shows an 
annulus of constant k. Notice that the power spectrum is not constant within 
this annulus. 


2.2 Redshift space and the wedge bias 

The real-space 21-cm signal—which is the typical output from 
reionization simulations—has no direction, and its power spectrum 
is isotropic. However, the real-space signal can never be observed. 
Any observation will see the signal in redshift space, where the ap¬ 
parent positions of emitters are distorted along the line-of-sight by 
their peculiar velocities (see e.g. Mao et al. (2012) for an overview 
of the effects of redshift-space distortions on the EoR 21-cm sig¬ 
nal). An emitter with a line-of-sight peculiar velocity tj|| and real- 
space position r will be translated to an apparent, redshift-space 
position s: 

1 + z ^ 

This distortion has two effects. First, it changes the contrast 
of the signal, making Ag (A) differ from A^{k). Second, since the 
translation happens along the line-of-sight, it makes the redshift- 
space signal anisotropic. The same holds for the redshift-space 
power spectrum, Ag (k); it is now a function of fi. In other words, 
to measure the spherically-averaged redshift-space power spectrum 
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Ag (fc), we need to average over an entire spherical shell in k space, 
including all values of jj.. Averaging only in the parts of a shell that 
are located in the EoR window will result in a biased measurement 
ofA?(fc). 

Figure 2, shows a 2D slice from a simulated 3D power spec¬ 
trum in redshift space. This figure makes apparent that the power 
spectrum is not isotropic and that averaging only outside the wedge 
will not accurately estimate Aa(fc). In the following sections we 
will investigate the magnitude of this error, which we will refer 
to as wedge bias, and define as [Ag „,i„dow(fc) - Ag (fc)]/Ag (fc), 
where is the power spectrum measured in the EoR 

window defined by some /imin- 



Figure 3. The mass averaged mean neutral fraction as a function of redshift 
in our simulations. 


3 SIMULATIONS 


To study the wedge bias on the 21-cm power spectrum, we use a 
set of semi-numerical reionization simulations carried out on top 
of a numerical A-body simulation. The A-body simulation was 
performed with CUBEP^M (Harnois-Deraps et al. 2013). This code 
calculates gravitational forces on a particle-particle basis for small 
distances and on a mesh for longer distances. We used 6912® par¬ 
ticles of mass 4 x 10^ Mq on a 13824® mesh, which was later 
downsampled to 600® cells. The total size of the simulation vol¬ 
ume was 500/h = 714 cMpc along each side. The minimum dark 
matter halo mass used in the simulation was 2.02 x 10® Mq. 

The reionization part of the simulations was carried out with a 
modified version of the semi-numerical code described in Choud- 
hury et al. (2009) and Majumdar et al. (2014). This code calculates 
the ionization state of the intergalactic medium by comparing the 
average number of ionizing photons entering in a cell with the av¬ 
erage number of neutral hydrogen atoms in it. A cell is considered 
ionized if it is possible to find a sphere of some radius (limited 
by the assumed mean free path of the ionizing photons) around it 
within which the number of ionizing photons exceeds the number 
of neutral hydrogen atoms. 

We use two different models for assigning ionizing fluxes to 
our dark matter halos. In the fiducial model, the ionizing flux is 
proportional to the halos mass, Mh'. 


iV.,(Mh) = iVion 


Mhflb 

rripllni 


(4) 


where Nion is the number of photons entering the IGM per baryon 
in collapsed objects and nip is the mass of a hydrogen atom. The to¬ 
tal number of ionizing photons is not conserved in this scheme due 
to the overlapping of ionized regions (Zahn et al. 2007). We tune 
the value of Nion at different redshifts to compensate for this and 
also to make sure that the resulting reionization history follows the 
same trend as the evolution of the mass averaged collapsed fraction 
with redshift. Our results (Section 4) are sensitive mainly to neutral 
fraction rather than redshift, and so we will present our results as a 
function of xhi for the remainder of the paper. 

The second model, which we will refer to as the “massive 
sources” model, has oc M^, with the proportionality constant 
tuned to give the same reionization history as the fiducial model. 
This model assigns higher fluxes to more massive sources, resulting 
in fewer and larger ionized bubbles. These two models were chosen 
to provide two very different examples of reionization topologies. 
The resulting reionization history (which is the same for both mod¬ 
els) is shown in Figure 3. 

After running the reionization simulations, we combine the 
density fields from the A-body simulations with the ionization 


Bmin —0.5 Bmin —0.95 



Xhi 

Figure 4. The difference between Ag(fc), A^{k) and ^j^j^^(A:). 
Dotted lines show the real-space power spectrum, A^{k), while solid 
lines show the redshift-space power spectrum Ag{k). Dashed lines show 
A®, window (fc) for Mmin = 0.5 (left panel) and ^^in = 0.95 (right panel). 
All curves aie for the fiducial model as a function of global neutral fraction. 


fields to get the 21-cm brightness temperature, making the simpli¬ 
fying assumption that the spin temperature is much higher than the 
temperature of the cosmic microwave background (CMB). Finally, 
we convert the output from real space to redshift space, using the 
methodology described in Jensen et al. (2013). All simulations use 
the cosmological parameters from WMAP five year data release: 
h = 0.7, flm = 0.27, Da = 0.73, = 0.0226 (Komatsu 

et al. 2009). 


4 RESULTS 

To illustrate the bias that occurs when measuring the spherically- 
averaged redshift-space 21-cm power spectrum in the EoR window, 
we calculate the power spectrum from our simulated 21-cm vol¬ 
umes, both for the full volume and for a window defined by some 
value of /imin. 

Figure 4 shows the power spectra for the fiducial model, in 
real space and redshift space both for an EoR window and for a 
full range of fj, values. It is clear from this figure that when we re¬ 
strict the measurement to certain ^ values, the wedge bias increases 
the difference between A® (fe) and A®(fc). Since the redshift-space 
distortion effect varies with redshift (see e.g. Mao et al. 2012), so 
will the wedge bias. The bias is most pronounced on large scales 
(small values of k) and, as expected, becomes more significant in 
the case of a smaller EoR window (higher fimin). Problematically, 
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M'min ~ 0-5 Prnin ~ 0.95 



Figure 5. The bias when measuring the spherically-averaged redshift-space 
power spectrum inside the EoR window, for a number of different k modes 
and two values of /imin- The results shown are for the fiducial model. 


14min — 0.5 P-min — 0.95 



Figure 6. Same as Figure 5, but for the massive sources model. 

increased bias at high redshift obscures the characteristic rise and 
fall signature in the 21-cm power spectrum. 

Figures 5 and 6 show the bias, i.e. [Ag window (fe) - 
Ag (A:)]/Ag (fe), when estimating the redshift-space power spec¬ 
trum in the EoR window, for our two different reionization models. 
The bias is positive early on, with the EoR window power spectrum 
over-estimating the redshift-space power spectrum by around 100 
per cent. At xm ~ 0.8, the bias turns negative and then becomes 
fairly insignificant in the later stages of reionization (less than 20 
per cent for all k modes even for ^min = 0.95). The change in sign 
is due to the increased anti-correlation between matter density and 
ionization fraction, which drives the change in the power spectrum 
anisotropy (see e.g. Mao et al. 2012; Jensen et al. 2013; Majumdar 
etal. 2013, 2015). 

The first, positive, peak in the bias is slightly more pronounced 
and occurs a little later in the fiducial model than in the massive 
sources model. The second, negative peak is almost absent in the 
massive sources model, especially on large scales. Overall, how¬ 
ever, the difference between the two models is not dramatic, despite 
the massive sources model being rather extreme. 


5 SUMMARY AND DISCUSSION 

When interpreting observations of the EoR 21-cm signal, one must 
always be wary of the difference between the fundamental real- 
space signal and the observable redshift-space signal. However, 
measurements using foreground avoidance, which sample only 


some values, also suffer from an additional bias since the power 
spectrum is anisotropic in redshift space. We have quantified this 
bias in two different reionization simulations for a range of red- 
shifts and spatial scales. We find that the bias is significant for 
higher redshifts and amplifies the difference between the real-space 
and redshift-space power spectra. Depending on how large a win¬ 
dow one can obtain, the redshift-space spherically-averaged power 
spectrum may be overestimated by around 100 per cent. In the later 
stages of reionization (xhi < 0.8), the bias instead results in an un¬ 
derestimation of the power spectrum. However, here the magnitude 
of the bias is smaller—less than 20 per cent on all k modes, even 
for /Tniin — 0.95. 

How can the wedge bias be accounted for in real-world obser¬ 
vations? When performing simulations to compare to observations, 
the situation is rather straightforward. As long as mock observa¬ 
tions from the simulations are produced with redshift space distor¬ 
tions included and using the same p cut as the observations, the 
wedge bias will be accounted for automatically. A more compli¬ 
cated situation is when comparing observations in the EoR win¬ 
dow to older simulation results from the literature. In this case, the 
wedge bias needs to be modeled. Our results indicate that the bias 
is not extremely sensitive to the reionization model, so it should 
be possible to model it to decent accuracy. In principle, it may be 
possible to determine directly if observed data falls in the low-bias 
regime by studying whether depends strongly on 

or not. However, this likely requires a large EoR window and low 
detector noise to be feasible (see Pober 2015 for a related discus¬ 
sion). 

Finally, we note that this study only presents a first, inexhaus- 
tive treatment of the wedge bias. We only show the effect for two 
different reionization models. The bias depends on the nature of 
the redshift-space distortions, which in turn depend on the reion¬ 
ization topology, and so the situation may be different for other 
reionization models. Furthermore, we have used the simplifying as¬ 
sumption that the intergalactic medium is heated at high redshifts, 
so that the spin temperature is much higher than the CMB tem¬ 
perature. This assumption may not hold at the highest redshifts, in 
which case the redshift-space power spectrum will look different 
(Ghara et al. 2015b). Since the wedge bias is highest at high red¬ 
shifts, it may be significantly affected by late heating. 

All the analysis in this paper was performed on data vol¬ 
umes with fixed redshift (“coeval volumes”). If the evolution of the 
signal—the lightcone effect—is taken into account, any measure¬ 
ment of the power spectrum will be done over a range of redshifts. 
The lightcone effect has been shown not to affect the anisotropy 
of the power spectrum at measurable k modes (Datta et al. 2014; 
Ghara et al. 2015a), and the wedge bias will be the same as for 
coeval volumes. 
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